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ABSTRACT 


i*/ri 

The variation of wind with altitude at 
low levels is hypothesized to be a universal 
velocity defect law; the exponential Sherlock 
formula fitted to observed data and extended 
to high altitudes is unbounded. Data from 
Hurricane Cleo 1964 and Dora 1964 is included 
with data from nine other storm happenings. 
Averaging of observed data as function is 
suggested in lieu of arithmetic averaging 

o 

of V; probabilities developed from V statistic 
(observed anemometer readings) need be trans- 
formed in probability to the linear V statistic 
to obtain high confidence levels at the extremes. 
Technical comments are solicited to promote 
math model development. 
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THE PROBLEM OF LOW LEVEL 


WIND DISTRIBUTION 

by J. H. Deese 
assisted, by 
Coralee C. Whisenant 

I. INTRODUCTION 

As a native Florida cracker, I can well remember the 
severe weather of exceptional thunderstorms and hurricanes. 

For instance, listening to the thump of our two story fra: e 
house as the wind lifted and dropped it again and again to 
its foundations or watching the shingles and framing of a 
section of the front porch fly apart and scatter with the 
sound of a freight train in a rumbling twister, or walking 
through the residue of a virgin pine forest where a hundred 
yard width swath of green timber had been twisted into match 
sticks. In those days of the early 20’ s, both timely weather 
warning and adequate structural design were non-existent. 
Today's Structural Engineering and architecture have as one 
of their prime purposes the design and construction of facil- 
ities to withstand the elements and to protect the people and 
equipments which they house. I am intensely interested in 
the development and usage of these practices. 

For the past fifteen years I have been intimately associat 
ed with the development of Service Structures at Cape Kennedy; 
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more recently, I am connected with the design and development 
of concepts employed at NASA Merritt Island Launch Area such 
as Arming Towers and the Vertical Assembly Building. 

In 1951/1952 the first "A frame” service structure was 
designed and built (Figure 1) . It was adapted from the exist- 
ing oil field portable drill masts. Wind loads were calculated 
on a constant horizontal wind from top to bottom; the structure 
height was 147 feet. Incidentally, this structure served 
Alan Shepherd, our first human astronaut. 

In 1956, the second ”A frame" structure 112 feet high was 
ordered (Figure 2) . In 1957 a third structure of the inverted 
U design, but still of the self raising leg type was purchased; 
its height was 134 feet (Figure 3). All of these structures 
were designed for reclining stowage in hurricane winds; moments 
from their erecting loads were greater than those imposed 
by hurricane winds in the guyed erect position. Overturning 
design moments without guying were calculated to occur at 
55 - 60 knots; this range had been selected as the maximum 
occurring afternoon thunder squalls. During the summer of 
1956 when Structure No. 1 was being used with a Redstone vehicle 
a 55 mile wind actually rocked the structure; (reminiscent of 
our old thumping homestead). On an expeditied basis, we 
proceeded to install a guying system to withstand 125 mph 
hurricanes $ our problem was ended. 


2 



When planning for Launch Complex 34 was initiated in 1958, 
design of Service Structure 34 was based on withstanding hurri- 


♦ 


^ cane force winds (125 mph) ; practice at that time introduced 

1 2 

the Sherlock exponential formula ’ . This structure was 310 
feet tall and weighed approximately 6,000,000 lbs (Figure 4). 
Wind effect was based on overall structure drag factor computed 
from wind tunnel tests; wind loads were computed from this area 
coefficient and the Sherlock wind profile with 1/7 exponent. 

Essentially the same process was repeated in the design 
of LC 37 Service Structure which was 303 feet high (top frame) 
and weighed approximately 10,000,000 lbs (Figure 5). Because 
of the open framing and cylindrical members, the effect of the 
Sherlock exponential profile was not significant as compared 
with the older practice of a constant value profile. This con- 
stant profile used the extreme recorded wind number. 

When feasibility studies and criteria development was 
initiated on MILA LC 39 (in June 1962), wind loads were found 
to seriously affect the Vertical Assembly Building structure; 
especially, those concepts which envisioned a completely 
enclosed structure offering weather protection to both the 
Saturn V space vehicle and the Launcher Umbilical Tower. A 
preliminary report^ comparing old and new practices presented 
» the problem; but the tight schedule to reach the Moon precluded 

extensive studies and attempts at extrapolation of weather 
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data from the 30 ft. U. S. Weather Bureau reference levels to 
higher altitudes. Our knowledge of the Nature of Low Level 
Winds (under 1000 ft. altitude) at Cape Kennedy was at that 
time almost nonexistent. 
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II. THE PROBLEM 


After a number of conferences with national authorities, 
the Architect Engineer^’ ® selected a velocity profile for hurri- 
cane exposure of the VAB at the MILA site (Figure 6) . This 
profile is representative of the best possible judgement that 
could be. made in the absence of comprehensive data. The factor 
of two and more by which the current ASCE practice differed 
from the older flat profile practice was most intriguing. Surely, 
if no failures existed in actual exposure history of the older 
designs, either the modern practice or the factors of safety as 
previously used are overly conservative. If one reviews the 
literature, }ie finds a lore of information on the characteristics 
of frames, members, and shapes when exposed to wind; but an 
almost complete dirth of information on the Nature of the Low 
Level Winds at specific locations such as Cape Kennedy and the 
Merritt Island Launch Area. 

As a comparison of the effect of various assumptions for 
expected extremal conditions, we plotted numbers from ASA and 

Q 

ASCE practices using reference velocities at 30 feet (Figure 7) ; 
the divergence of these at the higher elevations: prompted a 

O 

consideration of raising the reference level (Figure 8) . This 
alone tended to make the numbers easier to live with. The diver- 
gence in Figure 8 is still significant, but it is materially 
reduced from that of Figure 7. 
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Recalling occasions of previous storm passages, we recon- 
structed numbers according to memory of various old timers at 
Cape Kennedy and Patrick Air Force Base. Believing the exponen- 
tial formula to be the best expression, we compared these 
numbers as indicated on the curves of Figure 9. The exponent 
for the fit varies between .19 and .2; suprisingly, compares 
favorably with the .20 value discussed by Farber and Bell as 
found by Kamei. Of interest is application of a variety of 
exponents to the wind number 157 mph when the anemometer at the 
Jupiter lighthouse blew away in 1934 (?) . 

Some means was desired by which the variety of exponents 
applied to the Sherlock formula * could be compared (Figure 
10) o This we did by selecting a hypothetical cantilevered one 
foot width strip 500 feet high then computing moments and shears 
for each elevation from the ground up (Figure 11) . Of course 
they all converge at the top, but the relative factors as we 
get nearer the bottom are extremely divergent ; at most by a 
ratio of 2:1. Continuing the hypothesis as to a phantom can- 
tilevered mast that could withstand such loads, we computed 
the imaginary weights for such a steel structure (Figure 12) . 
The ratio of the numbers found is alarming; we ask ourselves 
"Are they truly representative?": - such as 1.34:1; 1„47:1; 
1,91:1; and 1.10:1. These precentages when applied to tonnage 
steel for large structures offers considerable promise that 
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savings on future construction can be effected and will be of 
no small value. The major requirement in such problem is that 
the Nature and Characteristics of Low Level Winds be better 
defined. Our curiosity had now felt the needle; my own personal 
desire for rationalization of the problem led into further inves- 
t igat ion. 



III. GENERALIZED DISCUSSION 


A review of flow theory^ ® ® (Figure 13) showed some 
promise. In viscous laminar flow both over flat plates and 
thru pipes, a maximum velocity is found in every instance. 

If the pipe be of infinite diameter, then like the plate, the 
maximum or free stream velocity exists in regions of the fluid 
where the viscous transported drag effect of the wall has become 
negligible. This condition is true for both laminar and tur- 
bulent streams having uniform average velocities. By analogy, 
there must be a region above the earth’s surface where the 
effects of its roughness on the wind becomes negligible. 

An early application of theory developed in the last 
century resulted in the Ekman spiral (Figure 14) ^ This 
function theoretically is limited by the gradient wind and the 
geostrophic wind. Various derivations by Lamb,^ Sutton,^ - ® and 
others make certain fundamental assumptions but the end form 

n 

of the equation is the same. That according to Lamb is given 
in Appendix I. 

If the maximum number exists, then an expression for the 
variation of velocity with height must exist such that 

K = U* K ' A 4r,j -ftt &*%*,*$] 

trtere: 

]//j ~ Xe/ae/Ay aAePeyaAa/? A j A?/ s/^e af rfsy/tPK 

V(yyp)“ X'X/j/ercepf dsj/MffeAe & f(X, ) *A 4 
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Such a function is bounded both by condition of ( V(pyp) — /)J 

in the free stream and by zero velocity at the surface. Although 
Sherlock recommends a constant velocity at altitudes above 1000 
feet and Kamei above 100 feet when using the exponential formula, 
the fact is that in 

V 


/ / \X r X depe/iJj on Socad/on, and 

K-K.( A) ■ 


/ 

as /j 


cx> 


and the function remains unbounded for x / 0 except at the 
surface. Therefore, it is an empirical expression of limited 
range. The problem remains of obtaining sufficient valid data 
to establish these limits. 

If one inspects flow turbulent theory as now practiced in 
wind tunnels, the Log Law of the Wall and the Velocity Defect 
Law are well confirmed. The Law of the Wall®’^ is exactly the 
Sherlock formula with 1/7 exponent; but its region of applica- 
tion is limited to the laminar sublayer. Outside this range, 
the Velocity Defect Law°> is receiving extensive investigation. 
Its fundamental expression is 



Definition of the point or range at which the transitions from 
Log Law to Defect Law occur has yet to be found. It is very 
likely that the limits of the transition functions can only be 
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So far these functions 


established by modern statistical theory, 
have eluded deterministic approaches. 
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IV. INSPECTION OF AVAILABLE DATA 


* 

For obvious reasons, such as knowing when to anchor down 
for wind extremes, Service Structures of LC 34 and LC 37 were 
instrumented with wind measuring devices. Sensing instruments 
are Bendix Friez Model 120 with frequency response of 60% per 
sec. Pen Recorders are of two types; Ester line Angus Model 
602 and Bendix Friez Model 141. Means is also provided for 
1” width MILAR tape recording for computer readout. Chart 
speeds of pen recorders are selective at 3 inches per hour and 
3 inches per minute. Installation of the aerovanes are shown 
in Figures 15 and 16. Lateral projecting instruments are 
mounted 8 feet from the structure (Figure 17) . Those on top 
of the structures are mounted at various heights (Figures 18 
and 19) above the structures. On top of each Blockhouse, 
instruments are also mounted (Figure 20) . 

On the pad in each complex, instruments are mounted on 
the Pad Light Pole (Figure 21) and on the Umbilical Towers 
(Figure 22). Unfortunately, there is some question as to the 

accuracy of these devices because of the structural influence. 

12 

Moses and Daubek reporting on a tower at Argonne National 
Laboratory, found that with 9 feet projection of the instruments 
serious errors occurred in the tower mounted Aerovanes due to 

A 

a tower influencing effect; both increases by 30% and decreases 
up to 50% were observed as well as direction deviations up to 
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40 degrees, Knothe and Callahan (Figure 23) have reported 

13 

tower influences of a similar nature. It appears that both 
wind tunnel model testing and consideration of aerovane prox- 
imity to the structure will be necessary to provide a means 
for correlation confidence in the observed values. Notwith- 
standing these faults, we have used the instrumentation plans 
of Figures 24 and 25 for recording eleven instances of extreme 
weather winds (Figures 26 thru 36) , The extreme maximum and 
minimum are plotted at each level as well as the average for 
that hour having the highest average and/or for the hour having 
the highest peak value. In either or both cases, the extremes 
within the period of record are shown. The data is presented 
in the following: 

Figure 26 Extreme Weather Winds, LC No. 34, 8 February 1964, 
0020 to 0120 hours 

Figure 27 Extreme Weather Winds, LC No. 37, 8 February 1964, 
0030 to 0130 hours 

Figure 28 Extreme Weather Winds, LC No. 37, 18 February 1964, 
1400 to 1500 hours 

Figure 29 Extreme Weather Winds, LC No. 34, 25 February 1964 , 
1500 to 1600 hours 

Figure 30 Extreme Weather Winds, LC No. 37, 25 February 1964, 
1420 to 1520 hours 

Figure 31 Extreme Weather Winds, LC No. 37, 28 April 1964, 

1330 to 1430 hours 
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Figure 32 Extreme Weather Winds, LC No, 37, 19 May 1964 to 
28 May 1964 

Figure 33 Extreme Weather Winds, LC No. 37, 27 August 1964, 

1545 to 1645 hours, Hurricane "CLEO" 

Figure 34 Extreme Weather Winds, LC No. 34, 9 September 1964, 
1400 to 1500 hours, Hurricane "DORA” 

Figure 35 Extreme Weather Winds, LC No. 37, 9 September 1964, 
1310 to 1410 hours, Hurricane "DORA" 

Figure 36 Extreme Weather Winds, LC No. 37 ^ 5 October 1964, 

1200 to 1300 hours 

These data include CLEO 1964 and DORA 1964. Trajectories for 
these are shown in Figures 37 and 38. Weather predictions of 
record excluding the hurricanes but regarding the other storm 
passages are shown in Table I, 

During Hurricane DORA, a curious buffetting caused baro- 
graph pumping; the barograph was located in a Southwest corner, 
second floor room of the E&O Building at Cape Kennedy (Figures 
39, 40, and 41). I had noticed a surging of the wind in 
Melbourne during the same time interval and on discussing it 
with the meteorologist he showed me the record (Figure 42) 
of the instance. The pen record shows quite clearly the jump- 
ing pressure; the solid portion was a continuous pumping. 

This buffetting at Cape Kennedy could have been caused by 
vortex shedding from the surrounding buildings^ since the 
winds were from the Southwest at that time. However, this 
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HAPPENING FORECASTS OF SEVERE WEATHER 
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does not explain my observation of the pulsing condition some 

30 miles away; both sound and eardrum pressure were evidence 

of the condition. I believe this to have been the nature of 

the hurricane, not a local condition. 

One of the most discussed storms subjected to analysis 

was Hurricane CARLA 1961 passage through Dallas, Texas, area 

(Figure 43). ^ The 1400 foot TV tower at Cedar Hill used by 

station KRLD had been instrumented for acceptance testing. 

This tower is located on the slope of a Knoll in an area of 

rolling hills (Figure 44). The entire period was reported by 

15 

Gerhardt et al. From his figures, the absolute magnitude 
of the horizontal wind (280 sec. average) was plotted against 
height. Hourly intervals are shown in Figure 45 together with 
their time ensembles; ten minute intervals are shown in Figure 
46 for that hour having highest winds. The curves indicate 
both stable flow and transient turbulence in the same storm. 

There is some question as to the validity of the instru- 
ment measurements at the 1400 foot level. Discarding these 
does not affect the form of the velocity function of height. 

In almost all instances, there is a vertical asymptote to V(h)* 
This data confirms that in this instance our model holds where 

^ = ^HYP jj ~f(h) Kyp; Q fit / A 7^] 

cm 4 l/tfyp ~ J 
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Again though, the influence of the structure on the observed 
velocity numbers is not known. The form of the function 
appears to be well defined. Also^ additional data is needed 
of the actual time *race rather than the 280 sec. averaged 
values. 
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V. ATTEMPTED ANALYSES 


For the period of September 20, 1963 thru October 1, 1963, 

1 fi 

an unusual wind condition existed. Point sets of recorded 
gusts occurring at each level are recorded as gust factor versus 
five minute averaged speeds in Figures 47, 48, and 49. Both 
the extreme and the mean show definite convergence with increas- 
ing velocity. The lowest level shows the greatest spread of 
observed gust factors, particularly at the lower speeds. The 
higher levels show a much lower spread of observed gusts at 
lower velocities (under 4 m/sec.). Cumulative distribution and 
density functions are needed for further analysis. 

For the 1964 happenings at LC 34 and LC 37, an aggregate 
of all recorded values was plotted. These aggregates are 
shown in Figure 50. On the plotted set of the left figure, all 
sets from the same data happening are connected; on the right 
figure the points are countably discrete. The sets are normal- 
ized to the average velocity squared value at each level. It 
can be seen that both the sets having values greater than 1 
and the set having values less than 1 indicate a convergence 
near 800 ft. height. Development of this could confirm a gen- 
eralization of approach suggested by the KRLD tower data. 

The cumulative distribution function of all aggregates 
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is plotted in Figure 51. The discrete density function of the 
observed values of the variable are shown below the plots. 


The Cumulative Function 

A- 




for values below 1.0 is accumulated in the negative direction; 
for values above 1.0, cumulation is in the positive direction 
on the abscissa number line. In Figures 52 and 53 the Cumula- 
tive Functions 

l 'V/u/c,) 

are similarly plotted. 

At this reading, the transform to continuous variables 
representative of these cumulative and density functions has 
not been completed. However, some conjectures can be made for 
these observations, such as: 

a. The extreme minimum of 




79. 81 //3; 8 3 i, 


j/ 7 / jsz y sxr 


for all cases of all levels is never greater than 0.8; and is 
bounded at 0. 

b. The extreme maximum of _ . 

for all cases of all levels is never less than 1.02; and is 
unbounded but is frequency limited with increasing values of 


' v/ivt; 


>/.o 
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Expectancy of any values greater than 12.0 will be very small. 

With collection of data from storms yet to come, more 
definitive analysis on these data and a further inspection of 
the extremes of point set functions of Figures 47, 48, and 49 
correlations can be obtained which will relate the density and 
cumulative functions of the aggregate to the velocities to be 
expected at each level. The completion of this effort will, 

I believe, be a beginning for the comprehensive analysis of 
Low Level Winds at Cape Kennedy and MILA. 
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VI. THE NEED FOR NEW APPROACHES 


14 

The current practices prescribed by ASCE paper No. 3269 
developed from the exponential Sherlock formula provide methods 
for estimating the winds and gust factors at altitude h; these 
average velocity and gust functions are limited to empirical 
ranges established by observed data. Doubt frequently exists 
as to the applicable limits of these ranges. Either reinspec- 
tion should be made of all available data for determining the 
validity of the Sherlock wind profiles or a new mathematical 
wind model should be developed along the lines already discussed. 
The deviation of the average value of the wind at any height 
from the asymptote to the velocity/height function at that 
height would then be representative of the maximum gusting that 
could occur for that particular storm. 

The frequency of gusting to any value within this range 
would be represented by the probability density function for 
that particular height. 

Because the model proposed herein would be developed from 
the extreme bounds, a correlation is needed between observed 
velocities for each class of wind storm that can exist and 
those forecasts of the meteorologist. The structural engineer 
is actually more interested in the extremes and their behaviour 
than he is in the predicted average value of velocity associat- 
ed with the happening. Records showing extremal predicted 
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values, classifications of the storm genesis, and the general 
synoptic conditions preceding and during the storm need to be 
developed. A new statistic is implied. 

The cumulative and density distributions of Figures 51, 

52, and 53 as well as their aggregates of Figure 50 are all 
plotted against height (h) and normalized V 2 . Most investigators 
use normalized V on a linear scale distributed variable. One 
may ask why we have not also used the linear relationship. The 
answer lies in elementals of mechanics and probability theory. 

The Aerovanes are momentum reaction machines; their signal 
output in most instances is linearly proportional to their 
speed 


SIGNAL = k N; N = rotational speed 
As rotating machines, they are not frictionless; the energy to 
overcome the friction is provided by the sensed wind. Their 
developed response characteristic (Appendix II) shows relation- 

“ “ / • v»- * ■ 

Therefore, the signal varies directly as the velocity squared. 

Since dynamic pressure is also proportional to product of 
density times V Q 2 , the signal output must vary directly 


as 


A 




Then 




or simpler, any observation with an aerovane or anemometer is 
2 

a V Q function. 
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In probability theory, the probability functions of 
variables of order 1 may not be used for their functional 
equivalents of order 'jjl, j (n ^ 1 )J without a transform in 
probability; also such transform must maintain a 1:1 point 
relationship for the countable or countably infinite sets. 
Simply, if a probability function is known for x, and if 

x = y 2 

then the p(x) must be transformed by mathematical operations 


to p(y) 


For the relationship 

2 - 




'0 0 

a transform method in probability for the two variables is 
developed in Appendix II. The illustrated case is that 
when p(N) is gaussian and the developed p ( C 0 V Q ) is Chi- 
squared. 

Also, in error analysis, the analysis of the primary 
variable must be limited to the observed quantities. There- 
fore, the governing distributions for probability develop- 
ment must be limited to deviates in the observation. 


Another example would be the average 





By the central limit theorem, Y and jTy^" are convergent 

with an infinite set but V never equals^V j a vg* i- n probablistic 

behavior. 
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Therefore, the natural statistic should be developed as 
f (V 2 ) . 

The resulting effect is, that for one function, the values 
of its density at the tail extremes differ by large orders of 
magnitude from that of its transform functions. Where extremal 
limits are of such importance as in the wind behavior at low 
levels, the expected values of extremes of particular happenings 
can be seriously misestimated, even with the most rigorous 
arithmetic . 

Surely, a new approach and development of new models of 
the low level wind behavior is indicated. 

The commentaries of all reviewers are welcomed to promote 
the search and development of these models. 
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Fig. 17 
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L.C. No. 37, RECORDS 19 MAY TO 28 MAY 1964 
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Looking northeast over the industrial area. 
FIGURE. 40 



FIGURE. 53 
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Looking east over the industrial area. 
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Fig. 45 
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FIGURE. 51 
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APPENDIX II 


WIND MEASURING REACTION MACHINES 
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